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Abstract 

The rigid pendulum, both as a classical and as a quantum problem, is an interesting system as 
it has the exactly soluble harmonic oscillator and the rigid rotor systems as limiting cases in the 
low- and high-energy limits respectively. The energy variation of the classical periodicity (r) is 
also dramatic, having the special limiting case of r —* oo at the 'top' of the classical motion (i.e. 
the separatrix.) We study the time-dependence of the quantum pendulum problem, focusing on 
the behavior of both the (approximate) classical periodicity and especially the quantum revival 
and superrevival times, as encoded in the energy eigenvalue spectrum of the system. We provide 
approximate expressions for the energy eigenvalues in both the small and large quantum number 
limits, up to 4th order in perturbation theory, comparing these to existing handbook expansions 
for the characteristic values of the related Mathieu equation, obtained by other methods. We then 
use these approximations to probe the classical periodicity, as well as to extract information on the 
quantum revival and superrevival times. We find that while both the classical and quantum peri- 
odicities increase monotonically as one approaches the 'top' in energy, from either above or below, 
the revival times decrease from their low- and high-energy values until very near the separatrix 
where they increase to a large, but finite value. 

PACS numbers: 03.65.Ge, 03.65.Sq 
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I. INTRODUCTION 



The problem of a point particle, restricted to a circular radius, and subject to a uniform 
gravitational (or electric) field, namely a rigid pendulum, is one of the most familiar of all 
problems in classical mechanics. The corresponding quantum version was first considered 
early in the history of quantum mechanics by Condon [jland is sometimes discussed in the 
more physical context of hindered internal rotations Q] -p in chemistry. The problem is also 
a staple of the pedagogical literature - 0] and is often seen as an example of perturbation 



3aagogica 
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theory methods |9J - especially in collections of problems in quantum mechanics 
- where is often described as the plane rotator in an electric field. Most treatments 
focus attention on aspects of the structure of the energy eigenvalues spectrum, with only a 
handful (see, e.g., [3]) examining the question of time-dependence of wave packet solutions. 

More recently, the realization that important information on longer-term quantum corre- 
lations, in the form of wave packet revivals and superrevivals, is also encoded in the energy 
eigenvalue spectrum, has focused new attention on the time dependence of physical and 
model systems which exhibit such interesting dynamics. In many such systems, initially 
localized states which have a short-term, quasi-classical time evolution, can spread signifi- 
cantly over several orbits, only to reform later in the form of a quantum revival in which 
the spreading reverses itself, the wave packet relocalizes, and the semi-classical periodicity is 
once again evident. Such revival phenomena have been observed in a wide variety of physical 
systems, esp ecially in Rydberg atoms [3], [3], and calculations exist for many other model 
systems [20]. 

The archetype of a one-dimensional model system for quantum revivals is the infinite well 
(where such revivals are exact) and a number of analyses 21] - j^j have provided insight into 
both the short-term and long-term behavior of wave packets. The rigid rotor, which like the 
infinite well, also has an n 2 quantum number dependence in the energy spectrum, behaves 
very similarly j^j], [2^, with neither system exhibiting higher-order correlations (such as 
superrevivals) due to the simplicity of their energy level structure. Harmonic oscillator 
systems, with a simple spectrum proportional to n, do not exhibit quantum revivals (or 
rather have T rev — > 00), or super-revivals, but Morse- type generalizations 29], with 

anharmonic terms present, have predictable revival structures, as confirmed in observations 
of molecular systems One-dimensional power-law potentials of the form V(k)(x) = 
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Vo\x/L\ k also provide examples |3l| of interesting quantum revival time behavior as a model 
parameter is varied, as they interpolate between the cases of the oscillator (k = 2) and 
infinite well (k — > oo) limits. 

For many such systems, which depend on a single quantum number, one typically ex- 
pands the energy eigenvalues (assuming integral values) about the central value used in the 
construction of a wave packet via 

E(n) « E(n ) + E'(n )(n - n ) + ~£>o)(n - n ) 2 + ~S'"(n )(n - n ) 3 + ■■■ (1) 

2 o 

in terms of which the classical period, revival, and superrevival times are given respectively 
by 

2nh 2ixh _ 2nh 

T= d ~W^o)\ ' \E"(n )\/2 ' super ~ \E"'(n )\/6 ( ) 

so that information on the long-term quantum correlations present in many wave packet 
solutions is encoded in the energy spectrum, just as is the short-term quasi-classical time- 
dependence. In typical experimental realizations of such systems, one finds r « T rev « 
T 

- 1 super ■ 

In this report, we will focus on the structure of the quantum revival (and superrevival) 
time for the problem of the quantum pendulum, extending earlier analyses of the approx- 
imate energy spectrum for this problem, in both the low- and high-energy limits, as well 
as making use of numerically obtained values for the entire spectrum. We are interested in 
such questions as (i) what is the revival time, T rev for the (anharmonic) oscillator limit of 
the quantum pendulum and (ii), how do the quantum revival and superrevival times vary 
(from both the low-energy oscillator and high-energy rotor limits) as one approaches ener- 
gies near the separatrix value. To probe such questions, we begin, in Sec. |H]by defining the 
generic problem of the quantum pendulum, and review standard expressions for the energy 
dependence of the classical periodicity of the pendulum, r(E) versus E. We then discuss, 
in Sec. IIII[ the quantum version of this problem, reviewing its equivalence to the Mathieu 
equation of the mathematical physics literature, and making contact with 'handbook' results 
for its periodic solutions from standard references. In Sees. IIVI and we discuss perturba- 
tive solutions to the quantum pendulum problem, in the high energy (free rotor) and low 
energy (oscillator like) limits respectively, deriving results up to 4th order in perturbation 
theory. We use these results to discuss the classical periodicity, as well as the quantum 
revivals and superrevivals in these two limiting cases. Finally, we use these earlier results 



as benchmarks against which to compare the numerical evaluation of the energy eigenvalues 
for the pendulum over the entire spectrum, focusing on the behavior of T rev and T super near 
the 'top' of the pendulum. Some of the mathematical details of the high-energy limit are 
discussed in a WKB-approximation in Appendix [XJ while the details of the perturbative 
expansion in the low-energy, oscillator-like limit are given in Appendix [B] 



II. CLASSICAL BACKGROUND 



For defmiteness, we will consider a point particle of mass \x (a notation used to avoid 
confusion with familiar angular momentum quantum numbers) restricted to a circle of ra- 
dius I, which therefore has a rotational inertia I = /i/ 2 . It is subject to either a uniform 
gravitational (mg) or electric (QE) force in the vertical direction, giving rise to a potential 
energy function of the form 

{mql (quantum pendulum) 
y ! (3) 

QEl (Stark effect for rigid rotor) 

(One can also consider a point electric dipole, p, in a uniform electric field, in which case 
V =pE.) 

We can rewrite the classical conservation of energy connection 

f(f) 2 -V cos W = S (4) 

in the form 

de 



(U / V E + V cos(9) (5) 
or _ 

7 de 

; = dt 6 

2 y/E + V cos(6) 

which can, of course, be formally integrated in terms of elliptic integrals (or evaluated 
numerically) to obtain 9(t) solutions. 

For energies satisfying —Vq < E < +Vq, the motion will be periodic with classical turning 
points at ±6 = ± cos _1 (— E/Vq) and the classical period is given by 

T(E)=2jir , = 
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where the extra factor of 2 comes from including both the 'back' and the 'forth' motions. 
For low energies, in the limit of small oscillations, with —Vq H E « +Vq or O << 1, the 



motion is oscillator like with an energy independent period given by r = 2ir/ y/Va/I. In the 
E — > +Vq (or separatrix) limit, the classical period diverges logarithmically. 

For the case of unbounded motion, where E > +Vq, one period corresponds more simply 
to a single revolution, with the result 

17 r +T de 

' [ l]Z= \2j^ y/E + VocosiO) 
and in the high-energy (rotor-like) limit, when E = IQ 2 /2 » +Vq, we find that 

1 2ix 2tt 



as expected. The expressions in Eqn. (JJJ) and (JSJ) for t(E) versus E can then be compared 
to the corresponding quantum values of r n (E n ) versus E n in the appropriate limiting cases. 

III. SCHRODINGER EQUATION FOR THE QUANTUM PENDULUM AND 
MATHIEU EQUATION 

The Schrodinger equation for the quantum pendulum can be written in terms of the 
moment of inertia, /, in the form 



&m + v(9)m = -^^P - v Q co S (e)m = Em . m 



The conventional change of variables, 9 = 2z, can be used to rewrite this as 



d 2 ij(z) 



+ 



ME 8IV ' 
H r^— cos(2z) 



h 2 h 2 



iP(z) (11) 



dz 2 

which is recognizable as a familiar equation of mathematical physics, 

^^ + [a + 2gcos(^)]^(^) (12) 

namely Mathieu's equation 3^|, j2] (with q — > —q in this case) and where we identify 

- 4IVo a - 8IE 

Because of the intrinsic parity of the potential, the solutions can be characterized as being 
even, ce r (z, q) or cosine-like, or odd, se r (z, q) or sine-like, for integral values of r, with r > 



(r > 1) for even (odd) cases. For the limiting case of q = 0, the solutions (in standard 
normalization 3^|, j^]) have the forms 



ce (z,0) = 1/V2 r = (14) 
ce r (;z,0) = cos(rz) r>l (15) 
se r (2i, 0) = sin(r2;) r>l. (16) 

For a fixed value of q, only for specific (or in more standard notation, characteristic) 
values of the parameter a will the solutions be periodic, with periods of ir or 2tt in the 
variable z, and these are denoted by a r ,b r for the even and odd solutions respectively. For 
each value of q, there are countably infinite number of solutions, labeled by r. For the 
physical case considered here, where the system must exhibit 2ir periodicity in the variable 
9, one chooses the even integer (r = 2m) ce 2m (<2,g) and se 2m (z,g) solutions which have 
period ir in the variable z. Standard mathematical packages (such as Mathematical) can 
quickly evaluate the characteristic values for a given value of q and we plot in Fig. 1 the 
pattern of characteristic values a 2m (solid curves for the even solutions) and 6 2m (dashed 
curves for the odd solutions) which are proportional to the allowed quantized energies, E, 
for different values of q, corresponding to increasing values of Vo- For E » +Vq (a >> q), 
the solutions are doubly degenerate corresponding to the free rotor limit: in this case, the 
real ce 2m and se 2m solutions are analogous to linear combinations of clockwise (exp(im9)) 
and counterclockwise (exp(— im6)) rotational eigenstates of the plane rotator. For — Vo < 
E « +Vo, the spectrum splits into equally spaced oscillator like states corresponding to 
the quantized version of small angle isochronous oscillations. 

Standard mathematical references (3^ , provide approximations for the characteristic 
values, Qjy ■ by ; CIS 8b function of q in both the q << 1 and q >> 1 limits, typically derived by 
expanding the solutions in appropriate Fourier-like series expansions. For the 'high energy' 
or rotor limit, with q << 1, one finds (for r > 7) that the approximately degenerate 

with 

a A=r> + ^^ + (5r2 + 7)g4 +■■■ . (17) 
br f 2(r 2 - 1) 32(r 2 - l) 3 (r 2 - 4) 1 ' 

and the difference between the characteristic values for even and odd solutions is known to 
satisfy 

a r - b r = 0(q r /r r ~ l ) as r -> 00 . (18) 



In the other limiting case, when q » 1 and the spectrum is approximately oscillator 
like, and one finds 

(p 2 + 1) (p 3 + 3g) (5p 4 + 34p 2 + 9) (33p 5 + 410p 3 + 405p) 
2g+2jV<? — gs ^ ^ ^y~2 ' ' ' ^ 

where p = 2ra + 1 which has the clear (n+l/2)foj oscillator dependence in low order. We will 
explicitly show, in Sec. and Appendix El how this expansion arises using perturbation 
theory, at least up to 4th-order, including all of the terms shown above. We note that 
handbook expansions exist to sixth-order (in both the low- and high-energy limits), so that 
one could improve on the explicit quantum mechanical perturbative calculations discussed 
below by using these results. 

While we will focus on approximate perturbative results in the next two sections, when 
we do compare to purely numerical results, we will make use of the following nominal set of 
parameter values, 

h = 2/i = I = 1 and V = 80 (20) 

corresponding to q — 160. This choice can be seen, from Fig. 1, to exhibit a large number 
of oscillator like states below the +Vq threshold. For these parameter values, Mathematical 
returns characteristic values with the equivalent of double-precision accuracy. 

IV. HIGH ENERGY ROTOR LIMIT 

In the limit of vanishing gravitational or electric field, the quantum pendulum reduces to 
a free rotor with Schrodinger equation and solutions given by 

-27%^ = ^° VmW with ^ )= dr iw * and E -=^f (21) 



with m = 0, ±1, ±2, .... The solutions with +m, — m for |m| ^ are degenerate, correspond- 
ing to the equivalence of clockwise and counter clockwise motions. Linear combinations of 
these solutions can be used to obtain eigenf unctions of definite parity as 

1 /y/2jr for m = 

cos(m^)/y / 7T for m > • (22) 
sin(m#)/ v /7r for m > 
The classical period associated with this zeroth-order result requires 



e m (0) 



dE m _ , h 2 m 2E m 

-d^ = E - = — = l— (23) 



so that 



2nh 



r«? es ff- = 2n/^/I (24) 



which is consistent with the purely classical result of Eqn. The corresponding revival 
time is 

2nh 4vrJ 

CV ~ \E'J2\ ~ h ■ {2b) 
These simple results can therefore be used to obtain zeroth- order approximations for the 
quantized energies, classical periods, and revival times in the large-energy rotor limit. We 
note that in leading order, the superrevival and all higher-order correlation times diverge 
due to the simple quadratic dependence on quantum number, just as for the one-dimensional 
infinite well. 

Using the standard set of values in Eqn. (|2()j) . we first exhibit the zeroth order energy 
eigenvalue predictions (compared to the 'numerically exact' results) in Fig. 2, and then 
compare the ri°' versus Em values to the classical periodicity in Fig. 3 in this same limit. 
Finally, in Fig. 4, we illustrate the predictions for the revival time, T rev versus E m . 

The effect of a constant external field (Stark effect) on the spectrum of the rotor is 

q n 

a well-known example of second-order degenerate perturbation theory 9] - |17] and the 
general result for the second-order energy shift is 

V) " h 2 {Am 2 _ 1) ^ 



for both even and odd states. Note that the first-order, and all odd-order shifts vanish for 
symmetry reasons. The special case of \m\ = 1 must be treated separately with the result 
that the +1,-1 states are split to this order jj| - jit]], [lfj], [l7|. More generally, it can be 
argued |6] that the degeneracy in a r ,b r is split in r-th order in perturbation theory, as in 
Eqn. (|18|) . When this result (Eqn. (|26j0 is written in terms of the parameters a and q, it 
reproduces the next-to- leading order of the 'handbook' expansion in Eqn. (|17|) to 0(q 2 ). 

This next-to-leading term then provides important corrections not only to the energy 
eigenvalues (especially for lower energy, smaller m states), but also to the quantum expres- 
sions for the classical period (via \E' m \) and the revival time (via |-E^J). Including these 
next-to-leading order results gives the second order predictions (shown as boxes) for these 
quantities in Figs. 2, 3, and 4. We note that the quantum revival times decrease from their 
high-energy value in Eqn. (J25|) as one approaches the critical E ~ +Vq value from above; 



9 



this is in contrast to the classical periodicity where the quantum r m values do seem to more 
closely follow the classical prediction of Eqn. (jHj) and exhibit signs of the expected diver- 
gence as E — > +Vq. This term also provides the first estimate of the superrevival time as 
one no longer has E'" = 0; this correction is somewhat akin to the introduction of a finite 



superrevival time in the finite well system 



3, 



35j | when the well depth is not taken to 



infinity, yielding a non-quadratic dependence on quantum number. 

While we are not aware of any corresponding explicit 4th-order calculation of the energy 
eigenvalues for the 'rotor in a field' system, one can make use of the perturbative expansion 
obtained by other methods in Eqn. (j!7j) to write the equivalent of the 4th-order perturbation 
theory result in the form 



£i 4) 



Vnl 2 







20m 2 + 7 



(27) 



k 6 L(4m 2 - l) 3 (4m 2 -4)_ 

and including these additional corrections does yield improved agreement with both the 
numerically evaluated energies (in Fig. 2) and in the prediction for the energy dependence 
of the classical period (in Fig. 3): it also sharpens the apparent 'dropoff' of the revival time 
as the energy approaches the separatrix value (as E —>■ +Vq from above) as illustrated in 
Fig. 4. Using the comparisons in Fig. 2 and 3, one can then better judge to what extent the 
predicted values for the revival times will likely reproduce the actual values for the complete 
quantum pendulum problem, as one approaches Vq from above. Henceforth, we will use the 
4th-order predictions for the revival (and superrevival) times in comparisons to numerically 
obtained values for the pendulum system. We note that the leading (in m) terms of the 
perturbation results in Eqns. ()21|). ()26|). and ()27|). as well as higher order ones in the detailed 
'handbook' expansion in Eqn. (fTTj) can be obtained in a very straightforward manner using 
the WKB approximation and we outline this procedure in Appendix [X] 

As a check on the standard assumption on the relative magnitudes of the various time 
scales in such problems, we can evaluate the lowest-order predictions for r, T rev , and T super , 
using the relations in Eqn. (j2J). The zeroth-order energy in Eqn. ()21|) gives the leading 
contributions to the first two quantities, while the second-order term in Eqn. (j2T?j) gives the 
leading term for T super and we find the hierarchy 

T supe r : T rev : r = (4m 3 /qf :2m: 1 (28) 

in the large m limit. We also recall that the next-to- leading corrections for each of these 
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quantities imply that r increases as E — > +Vq, while both T rev ,T super decrease so that the 
time scales are compressed near the separatrix. 



V. LOW ENERGY HARMONIC OSCILLATOR LIMIT 

In order to study the low-energy or oscillator-like limit of Eqn. (JTOJ), we write / = fil 2 
and expand the potential energy function for 9 « 1, giving 

h 2 d 2 tb(9) f 9 2 9* 9 6 \ , s 

Collecting only the constant and quadratic terms from the potential energy, and writing x 
01 (again, in the spirit of a small angle approximation), we can write the exact Hamiltonian 
in the form 



r=2 



and we will treat the various H2r terms in H' as perturbations. We next equate 



so that the zeroth order energy eigenvalues are trivially given by 

= ( n + l/2)Hu -V = (n + l/2)h v / Vo/I - V . (32) 

Expressed in terms of a and q, this reduces to 

a = -2q + 2(2n + l)^q (33) 

which are the two leading terms of the 'handbook' expansion in Eqn. (|19|). This harmonic 
oscillator approximation is compared to the numerically obtained quantum pendulum eigen- 
values in Fig. 2 (as crosses) as the leading (zeroth) order prediction, for E < +Vq. The 
classical periodicity from Eqn. (J2J) is given by 

2uh 2tt 



\E'\ uj 



(34) 



and we superimpose these zeroth order predictions for (E n , r n ) on the classical prediction of 
Eqn. (JZj) in Fig. 5 (as crosses). 
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The leading correction in this approach comes from treating H 4 using first-order pertur- 
bation theory to write 



E n ] = (n\H4\n) = -^(n\x 4 \n) = 



°" -X2n 2 + 2n + l) 



(35) 



using elementary textbook results for the expectation value of x A in oscillator states. As 
always, the evaluation of the expectation value of powers of x is most easily done using 
raising and lower operator formalism, with 



h 



x 



(A + i f ) (36) 



2/io; 
where 

A\n) = y/n\n - 1) and A ] \n) = Vn + l\n + 1) (37) 

and we use this method in all our calculations below. Using the identification in Eqn. (|31|h 
we find that this reduces to 

EV = -^ ! (2n* + 2n + l). (38) 

In terms of a, q, this contributes a Vo-independent correction of the form 

fl(1) = _ (^+4^ + 2) = (p* + 1) (3g) 
8 8 

where p = 2n + 1, which is the g- independent term in Eqn. 

This correction, arising from the anharmonicity of the potential, not only provides an 
improved prediction for the energy-period relationship in Fig. 5, but is also important as it 
gives the lowest-order prediction for the revival period in the low-energy limit, as in Refs. 2^ 
and 29]. The revival time is given by 

2ttH _ 327r/ 

rev ~ \Kj2\ ~ ~h~ ( } 

which is independent of Vq and 8 times the lowest-order result for the revival time for the 
pure rotor system from Eqn. ([25)1 . Thus, the revival time in the low-energy limit is predicted 
in terms of physical parameters of the systems, by the very specific anharmonicity arising 
from the form of the cos(6 l ) potential. 

We show this first non- vanishing prediction for the revival time in the low-energy, oscilla- 
tor like limit in Fig. 6 (as diamonds.) At this order, the superrevival time is formally infinite 
since the spectrum is still approximately quadratic in n. 
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At next order, one obtains two contributions, one from including the effect of Hq in first 
order perturbation theory, and a second from using H4 in second-order. The results can be 
written in the forms 

= (n\H B \n) = +^j E (n\x e \n) = ^ (J^ 5(4n 3 + 6n 2 + 8n + 3) (41) 



and 



(n\H 4 \j){j\H 4 \n) 

(0) 77,(0) X 



Vo \ 1 v- (w|x 4 |j)(j|x 4 |n) , . 

4!Z 4 J fruZf (n-j) [ } 



4!Z 4 



2 



1 



(34n 3 + 51n 2 + 59n + 21) 



Again, using to = \/Vq/I we then find the total second-order contribution to be 

E (2) = E (2,a) + bM = g /T (2n 3 + 3n' + 3n + 1) 



P V 512 < 43 » 

which agrees (when expressed in terms of a and q) with the expansion in Eqn. ()19|) and 
improves agreement with the 'exact' (numerically evaluated) pendulum values as shown in 
Fig. 2. The (E n , r n ) values including these corrections are also plotted in Fig. 5 and similarly 
show improved agreement in the E < +Vq regime. In analogy to the high-energy rotor case, 
the lowest-order prediction of a fixed value of the revival time in this limit, in Eqn. (|4()jl is 
changed in the next order so that T rev decreases as E — > +Vq from below. 

We can proceed in an identical (if increasingly more difficult) fashion to find the 3rd and 
4th order corrections, namely 

(3) h A (5n 4 + 10n 3 + 16n 2 + lln + 3) 

n ~ ~VqP 8192 ( } 

(4) h 5 (66n 5 + 165n 4 + 370n 3 + 390n 2 + 225n + 53) 

and some of the details of their calculation are presented in Appendix [BJ In each case, the 
results match the 'handbook' expansion from Eqn. ()19j1 : they also improve agreement with 
numerically obtained values (as in Fig. 2) and with the classical periodicity (in Fig. 5), as 
well as steepening the 'dropoff' of the predicted quantum revival time in this limit as shown 
in Fig. 6. We see that the quantum revival times decrease (at least initially) as one goes 
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up in energy, in contrast to the quantum periods (r„) which increase in accordance (again, 
initially) with the classical prediction as one approaches the classical divergence in r. 

We can also check the hierarchy of time scales, as done at the end of Sec. I1V| in the 
low-energy limit. In this case, the lowest-order predictions for r, T rev , and T super come from 
Eqns. and respectively and give the simple ratios 

T super : T rev : r = (8^) 2 : (8^) : 1 (46) 

and higher-order corrections lead to increasing (decreasing) values of r (T rev ,T super ) as E —* 
+Vq as in the high-energy rotor limit, with the various time scales being 'compressed' near 
the separatrix. 



VI. CLASSICAL PERIODICITY, REVIVAL AND SUPERREVIVAL TIMES FOR 
THE PENDULUM 

We have derived perturbative expressions for the quantized energies in both the low- and 
high-energy limits, and can therefore evaluate approximate expressions for all of the relevant 
time scales in Eqn. (j2J), as well as judging their relevant range of validity by comparing to 
'exact' values as in Fig. 1. In order to examine the behavior of the classical periods, and 
revival and superrevival times, near the separatrix (for E « +Vo) where neither approxi- 
mation scheme is applicable, we need to use numerically obtained values. With the energy 
spectrum, E n versus n, obtained in this way, we can approximate the required derivatives 
with respect to quantum number by making associations such as 



dE n AE n „ „ , d 2 E n A^E n ^ „ 

-r 1 -r^ = E n+1 - E n and —^^——^ = E n+1 -2E n + E n . 1 47 
an An an^ An z 

with similar expressions for higher order derivatives. However, because of the degeneracy 
between even and odd rotor-like states in the high energy limit, if we simply enumerate all 
of the allowed states in some standard ordering, we find vanishing 'derivatives' in the free 
rotor limit which are artifacts. In the construction of a general, high-energy wave packet, 
one would find approximately equal contributions from both the (almost exactly degenerate) 
even (ce2 m ) and odd (se2 m ) solutions which would be only 'counted' once. To avoid this 
difficulty, we can classify the states by their parity and evaluate the approximate derivatives 
in Eqn. (|4Tj) separately in both the even and odd cases. For the low-energy, oscillator limit, 
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this means that we 'sample' every other state and so, for example, obtain first-derivative 
type terms which are a factor of two too small (since effectively we have used An = 2 instead 
of An = 1), with a similar factor of 4 for the second derivative term. This procedure is not as 
artificial as it may seem, since one can imagine constructing wave packets of definite parity, 
even or odd, at appropriate points of symmetry for the pendulum (either at 6 = or 9 = 7r) 
which would then have contributions only from the even or odd parity eigensolutions. 

Following this procedure, we construct approximations for the first derivatives by gener- 
alizing Eqn. (|47|) slightly by associating the 'average' energy 

1 A E 

K=-(E n+1 + E n ) with — — ^ = E n+ i — E n = AE n (48) 

2 An 

so that we associate 

— - , _ 2irh . 

E n and r n = —— . (49) 



\AE. 



n | 



We plot, in Fig. 7, both the even-parity (crosses) and odd-parity (diamonds) values for 
(E n ,T n ) obtained in this way. (Note that the values for E < +Vq are shown on an appro- 
priately scaled (factor of 2) axis, compared to those values for E > +Vo.) We note that the 
'exact' data do closely follow the classical curves from Eqns. (JZJ) and (|8|). and do exhibit a 
large, but finite peak at E ~ +Vq, in both sets of data. The lack of a true divergence is 
understandable given the quantized nature of the energy eigenvalues and has been discussed 



which 



in terms of the explicit time evolution of wave packet solutions as described in Ref. 
also outlines simple uncertainty principle arguments. 

Turning now to the question of the behavior of the revival time as E — > +Vq (from either 
above or below), while there is no classical prediction comparable to that for r(E) versus 
E, we can make use of the 4th-order perturbative results to evaluate T rev = 2irh/\E"/2\ 
versus E n which will give good representations in both the low- and high-energy limits. 
Using the 4th-order results plotted in Figs. 3 and 5, we produce the limiting predictions for 
T rev versus E n , shown as solid curves, in Fig. 8. The numerically obtained 'exact' results, 
using the prescription in Eqn. ([47)1 . are also shown in Fig. 8, again with the low-energy 
'data' (E < +Vo) shown scaled, now by a factor of 4. Once again, the agreement with 
the perturbative predictions in the two limits is good, with T rev clearly decreasing from the 
appropriate E « +Vo and E » +Vq limits until very near E « +Vo where it reaches a 
large, but finite value. Similar results can be seen for the superrevival time, T super , comparing 
the predictions from 4th-order perturbation theory in both the low- and high-energy limits 
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with discretized derivatives of the numerically obtained results (appropriately scaled for 
E < +Vo), yielding plots which are quite similar to Fig. 8. 

VII. DISCUSSION AND CONCLUSIONS 

We have examined the quantized energy eigenvalue spectrum of the quantum pendulum 
problem, focusing on information about the classical periodicity and quantum revivals and 
superrevivals. Using 4th-order perturbative approximations for the low-energy (oscillator) 
and high-energy (rotor) limits of the pendulum, we have found good agreement with nu- 
merically obtained eigenvalues. This has allowed us to derive expressions for the revival 
and superrevival times in these limits, as well as to establish limits on their validity. The 
revival times for the pendulum decrease from constant values, proportional to h/I, as one 
approaches the separatrix from both above and below. Near the separatrix, we have used 
numerically obtained values to show that the classical period exhibits a large, but finite 
value at E « +Vo, while the revival times also increase (very near this value). 

This problem is an interesting example of a quantum system which exhibits a rich variety 
of classical periodicities and quantum revival and superrevival times, all of which are per- 
turbatively calculable in limiting well as classically divergent behavior at a special 

value, which is 'softened' by quantum effects. The behavior of the time scales as they near 
the critical value at the separatrix, where they can be clearly seen to compress towards 
each other as E — > +Vq, differs from the pattern seen in other systems with quantum re- 
vivals. One can speculate on possible experimental realizations of such behavior in systems 
exhibiting hindered internal rotations [J - Q]. 

The energy level structure of the quantum pendulum is, to a large extent, determined 
by the numerical (and dimensionless) value of q and we can discuss order-of-magnitudes 
estimates for q for some model systems, and how they scale with the physical parameters 
involved. In the classical limit of a macroscopic pendulum system under the influence of 
gravity, with m ~ 1 kg and I ~ 1 m, we find q ~ 10 70 so that there is the expected huge 
number of oscillatory states below the separatrix and one is clearly in the correspondence 
principle limit. At the other extreme, for an electron (m e ) restricted to a radius characteristic 
of that of a carbon nanotube (L ~ lnm), under the influence of gravity, we find that 
q ~ 10~ 18 and the system is automatically in the high-energy or free-rotor limit. In this 
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case, the revival time should be given to an excellent approximation by Eqn. (|25|). with 
no corrections due to the external potential, as well as having T super — > oo. For the more 
interesting case of an electron (or ion) subject to an electric field of order E = 100 V/m, we 
find q ~ lO~ 6 (M/m e )(E/EQ)(l/l nm) 3 so that for electrons in such fields the high-energy 
free rotor limit is also applicable. On the other hand, for ions (M/m e > 10 4 ) in strong 
fields (E > 10 5 V/m) in larger diameter (L > 10 nm) geometries, one has q ~ 10 4 and 
an interesting number of vibrational states could be present. This range is perhaps more 
typical of that explored in studies of hindered internal rotations Q| We have found 
that the revival times, in both the high-energy (Eqn. (|25jl ) and low-energy (Eqn. ()4()jl ) 
limits scale as T rev oc I/h, which gives numerical values in the range, T rev ~ (100 fs — 
1000 fs) (M/m e )(l/lnm) 2 , independent of the pendulum potential, V . 

A simple extension of this system would be to consider a particle confined to a circular 
(radius I) cylinder of length L, subject to a field perpendicular to the central axis. This 
separable system would then have quantized energies consisting of the quantum pendulum 
values discussed here plus a term given by E% = h 2 k 2 /2fiL 2 corresponding to the quantized 
motion along the cylinder. The pattern of classical periodicities and quantum revivals in 
systems characterized by two quantum numbers has been discussed j3] , {stJ . A more careful 
study of this case might be motivated by its possible realization in the geometry of carbon 
nanotubes. A related two-dimensional generalization would be to consider the point object 
restricted to the surface of a sphere which has a richer set of possible classical motions. 
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APPENDIX A: WKB APPROXIMATION SOLUTIONS TO THE MATHIEU 
EQUATION 



We briefly describe in this Appendix a simple WKB-based approximation which gives the 

leading order (in m) results for the (approximately degenerate) characteristic values, a r ,b r , 

in the high-energy limit (rotor-like) of the quantum pendulum problem, as an expansion in 

powers of Vq. In this context, it reproduces the leading m behavior of well-known 'handbook' 

results in Refs. jj^J and (^J for the characteristic values a m , b m in terms of the parameter 

q, in a very straightforward way 

We first rewrite the classical energy conservation connection 

2 

^L-V cos(9) = E (Al) 

in the form 

Pe = V2I^/E + V cos(6) . (A2) 
In the spirit of a WKB approach, we quantize the integral of po in the form 

/+7T P+TT 
+ V Q cos(fl) d6 = / p de = m7rh. (A3) 
-7T J — TV 

For the physical situation of the quantum pendulum, we would restrict ourselves to even 
values of m, but for comparison to the general result inspired by Mathieu equations, we 
allow any integral m values. In terms of the Mathieu equation parameters in Eqn. (jl3j) we 
can write this in the form 

2mvr = f + — cos(6)) d6 (A4) 



a 

One can expand the integrand as a series in q/a and perform the angular integrals, leaving 
only the even terms present, to obtain 

/ q 2 15g 4 105g 6 \ 

m= ^ a \ l -^-^-2w + ---) (A5) 

or 

2 ( q 2 15q 4 105g 6 \ 2 

a — m 1 + ■■■]. Ao 

V 4a 2 64a 2 256a 6 J K J 

This can be solved recursively in powers of q 2 to obtain the approximation solution 

q 2 5g 4 9g 6 
2^? + 32m 6 + 64m 10 



m i + jL + ^ + ^L i+ .... (A7) 



and this form agrees with the handbook results of Refs. 32] and 



3- 



One can then easily 



extend this result (again to leading order in m) to quite high order. 
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APPENDIX B: THIRD- AND FOURTH-ORDER PERTURBATION THEORY 
RESULTS IN THE OSCILLATOR LIMIT 



For the third-order correction to the low-energy, harmonic oscillator limit, we need to 
evaluate three distinct contributions, namely, (i) the H$ term in first order, (ii) the H$ and 
H4 terms (in both orders) in second order, and (iii) the H4, H4, H4 terms in third order. 
For example, the first contribution is 



/ 8 8! 



(n\x \n) = — 



h 2 



5(14n 4 + 28n 3 + 70n 2 + 5Qn + 21) 
8!2^ (B1) 



The expression for the third-order energy shift is presented in some elementary texts, but 
we reproduce it below for completeness, specifically 

F m _VV (n\H'\j){j\H'\l){l\H'\n) _ ~ v (n\H'\j)(j\H'\n) 

Z^Z^/ p (o) p(o)wp(o) p (ok \ n \ n \ n iz^ / p (o) p (o)x 2 y DZ > 



Including all three sets of terms, we obtain the result shown in Eqn. 

For the increasingly complex matrix algebra involved in the evaluation of such terms, we 
make use of symbolic manipulation programs (in our case Mathematical) and write 



(i\x\j) = J 2^<^ + i f |j) = y ^ + vT+l^ij+i] (B3) 

and evaluate the required matrix elements of higher powers by repeated matrix multiplica- 
tion. We note that in third-order (and beyond) in perturbation theory, there are consistently 
terms including factors of the form (n\H'\n). At third-order, and beyond, in intermediate 
steps of the calculation, the individual expressions one obtains may be of higher order (in 
n) than are present in the final result, and the presence of such cancellations can be used as 
a check on the internal consistency of the calculation. 

Finally, for the 4th-order terms, we require the contributions from (i) H\q in first order, 
(ii) H$, (both orderings) in second order, (iii) H 6 , H 6 , also in second order, (iv), H 6 , H4, H± 
(three orderings) in third order, and (v) H4, H4, H4, H4 in 4th-order. The explicit expression 
for the energy shift in 4th order, in terms of unperturbed wavefunctions and energies, is less 
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often seen (see, e.g., 
£l 4) 



so we also reproduce it here for completeness. 

(n\H'\j){j\H'\l)(l\H'\r)(r\H'\n) 



/ ; , 

j'^n Ij^n r^n \ 



(0) 
re 



£ 7 {0) )(^ 



(0) 



j^n i^n l"^ 



, ^)(4 0) - 

(n|^|j)(i|^|0(#'|n) 



E. 



{oh 



(0) 



^ 0) )(^re 0) 



(0)^2 



-(n|ff'|n>X;E 



j^n l+n ~ Ej ) {E n - E { 

(n\H'\j)(j\H'\n) 



(0) WP (0) 



+ (n\H'\n)(n\H'\n) 



(0)^ 



(B4) 



E 



(o) 



(n|^"'|r)(r|^"'|n) 



^ (E {0) -E { 

As a cross-check of our more automated approach to the 3rd and 4th order calculations, 
we have followed the same approach for the 1st- and 2nd-order results and confirm that we 
reproduce the values in Eqns. (|3*H|) and (|43|) . Perhaps more interestingly, we also note that 
there are three simple extensions of the 'basic' oscillator, defined by the zeroth-order poten- 
tial V(x) = fiu 2 x 2 /2, which can be used as benchmarks for such higher order calculations. 
First of all, the addition of a constant perturbing potential term, Vq, should only affect the 
energy spectrum through the first-order correction En = (n|V^|n) = V ' and one can easily 
check that all higher order corrections vanish in this case. 

The inclusion of a linear term of the form V'(x) = —Fx can be solved exactly by com- 
pleting the square giving 

2 m9 



(n\H'\j)(j\H'\n) 



E (E (o) _ 



1 OO 1 Q 

V(x) = -jjw x — Fx = -/jw 



X 



2^ 



V(x - x Q ) - Eq 



(B5) 

Fx perturbation 



to give E' n = E n — F 2 /2/j,/lv 2 = E n — Eq. Thus, the inclusion of an H' 

(2) 

should find a non- vanishing result only for En and we confirm this. 

Finally, the inclusion of a small quadratic correction to the original oscillator potential 
by V'{x) = Xx 2 leads to a simple redefinition of the oscillator frequency uo' 
and the result 



'I + 2\/hio 2 



E' 



[n 



l/2)hw' = E n il + 



2X 
ftjjj 2 



1/2 



(B6) 



which can be expanded as an ordinary power series in A. One can then compare this result 
to that obtained by evaluation of the first- through 4th-order perturbation results, including 
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those of Eqns. (|B2|) and (jB4|) and we find complete agreement. 
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FIG. 1: Plots of the characteristic values for the Mathieu equation, a2 m (for even solutions, solid 
curves) and 62m (for odd solutions, dashed curves) versus q for the quantum pendulum. The dotted 
lines correspond to a = ±2q or E = ±Vq. Note that for a given value of q (at least for a << q) 
that the gap between the lowest energy state (lowest solid curve) is roughly one-half of the spacing 
between solid and dashed curves, corresponding to the zero-point energy in the oscillator limit. 
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FIG. 2: Plots of the difference between the perturbative approximations for the energy eigenvalues 
and the 'exact' values (from direct numerical evaluation) in various orders of perturbation theory. 
The parameter values in Eqn. I|20jl are used. For the E > +Vq = 80 cases (using the high energy 
rotor limit) the zeroth, second, and 4th-order expressions in Eqns. (|21|). (|26|). and (|27jl are used. 
For the E < +Vq case, the low-energy oscillator approximations in Eqns. (J32J), ((581) . l(15|) . ((H)) . and 
P51) are used. 
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FIG. 3: Plots of t{E) versus E for the classical pendulum from Eqn. (JHJ) (solid curve) and of the 
equivalent quantum values, r m = 2i:h/\E' m \ versus E m using zeroth- (from Eqn. (|2"Tj) . as crosses), 
second- (from Eqn. as squares), and 4th-order (from Eqn. ((771) . as stars) perturbation theory. 
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FIG. 4: Plots of the revival time T rev versus E using T rev = 2i:h/\E'J n /2\ in zeroth-, second-, and 
4th-order perturbation theory (with the same notation as Fig. 3.) 
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FIG. 5: Plots of t{E) versus E for the classical pendulum from Eqn. (J7J) (solid curve) and of the 
equivalent quantum values, r m = 2irh/\E' m \ versus E m using zeroth (from Eqn. ()32|) . as crosses), 
first (from Eqn. (J38|). as diamonds), second (from Eqn. (|43j) . as squares), third (from Eqn. (|44j) . as 
bursts), and 4th order (from Eqn. I|27|). as stars) perturbation theory. 
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FIG. 6: Plots of the revival time T rev versus E using T rev = 2irh/\E'^ n /2\ in first-, second-, third-, 
and 4th-order perturbation theory (with the same notation as Fig. 5.) 
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FIG. 7: Plot of r n = 2nh/\ AE n \ versus E n for the parameters in Eqn. (|2Uj) for even (crosses) and 
odd (diamonds) energy eigenvalues. The classical expressions for r{E) versus E from Eqns. (J7|and 
are used. 
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FIG. 8: Plot of the revival times T rev = 2ttH/\ A^ 2 >E n /2\ versus E n for the parameters in Eqn. ()2U|) 
for even (crosses) and odd (diamonds) energy eigenvalues. The solid curves are the predictions 
using 4th-order perturbation theory from the respective limiting cases. 



27 



